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ABSTRACT 

(N| ; We present HST Faint Object Camera ultraviolet (A = 2300 A, 1400 A) images 

' of three ultraluminous infrared galaxies (ULIGs: L ir > 10 12 L©)[] selected from the 

■ IRAS Revised Bright Galaxy Sample. The purpose is to estimate spectral energy 

distributions (SEDs) to facilitate the identification of similar objects at high redshift 

ff) . in deep optical, infrared, and submillimeter surveys. All three galaxies, VII Zw 031 (= 

IRAS F05081+7936), IRAS F12112+0305, and IRAS F22491-1808, were well detected 

Q\ \ with the F220W filter. Two of the three were marginally detected with the F140W 

filter. The fluxes, together with ground-based optical and infrared photometry, are 

i-£h ■ used to compute SEDs over a wide wavelength range. The measured SEDs drop from 

the optical to the ultraviolet, but the magnitude of the drop ranges from a factor of ~ 3 
in IRAS F22491-1808 to a factor of ~ 100 in VII Zw 031. This is most likely due to 
different internal extinctions. Such an interpretation is also suggested by extrapolating 
to ultraviolet wavelengths the optical internal extinction measured in VII Zw 031. 

. £h ! ET-corrections are calculated to determine the colors of the sample galaxies as seen at 

high redshifts. Galaxies like VII Zw 031 have very low observed rest-frame UV fluxes 
which means that such galaxies at high redshift will be extremely red or even missing 
in optical surveys. On the other hand, galaxies like IRAS F121 12+0305 and IRAS 
F22491-1808, if seen at high redshift, would be sufficiently blue that they would not 
easily be distinguished from normal field galaxies, and therefore, identified as ULIGs. 
The implication is then that submillimeter surveys may be the only means of properly 
identifying the majority of ULIGs at high redshift. 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by AURA, Inc., under NASA Contract NAS 5-26555. 

2 Z/i r = L(8 — 1000 (im) computed from the flux in all four IRAS bands according to the prescription given in 
Perault (1987; see also Sanders & Mirabel 1996). Throughout this paper we assume that Ho = 75 km s _1 Mpc -1 
and that Qq = 1. 
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Subject headings: galaxies: formation - galaxies: fundamental parameters - galaxies: 
photometry - galaxies: starburst - infrared: galaxies 

1. INTRODUCTION 

Dissipation plays a major role during galaxy formation. It may be the most important 
process that determines the fundamental properties of the visible parts of galaxies (Blumenthal 
et al. 1984). Ultraluminous infrared galaxies (ULIGs) are the best candidates we know about for 
galaxies that are currently undergoing dissipative collapse (Kormendy & Sanders 1992). Their 
huge infrared luminosities (L- u > 10 12 L Q ) result from circumnuclear starbursts and/or active 
galactic nuclei (AGNs) surrounded by dust. These galaxies contain large amounts of molecular 
gas (e.g., Sanders et al. 1991), with a large fraction of this gas concentrated close to the nucleus 
- typically Mh 2 ~ 1 x 10 10 Mq in the central lkpc diameter region (e.g., Solomon et al. 1997; 
Scoville et al. 1997) - which presumably serves to fuel the nuclear activity. The morphologies of 
local ULIGs suggest that the collapse was triggered by a recent merger (Joseph & Wright 1985; 
Sanders et al. 1988a; Clements et al. 1996; Murphy et al. 1996). 

Kormendy & Sanders (1992) argue that ULIGs may be giant elliptical galaxies in formation: 
they have solved the well-known problem (Ostriker 1980; Carlberg 1986; Gunn 1987; Kormendy 
1989, 1990) of producing the observed high core densities of ellipticals via an ongoing dissipative 
collapse and starburst. Preliminary indications from the still-small number of ULIGs that have 
been observed at high resolution suggest that they are consistent with the core fundamental plane 
correlations (e.g., Kormendy & Djorgovski 1989) of elliptical galaxies (Doyon et al. 1994). It has 
also been argued that most quasars and other powerful AGNs began their lives as ULIGs (Sanders 
et al. 1988a,b; for a recent review of these arguments, see Sanders & Mirabel 1996). Therefore 
low-z ULIGs, although rare, may provide an excellent opportunity to study the processes by which 
ellipticals and AGNs formed. But before we can apply this hypothesis at the distances where most 
present-day ellipticals formed and where quasars were most abundant, we need to be able to find 
ULIGs at high redshifts. 

Ground-based followup observations of objects discovered by the Infrared Astronomical 
Satellite (IRAS), and more recently by the Infrared Space Observatory have been productive in 
identifying ULIGs. However, relatively few objects have been identified from these surveys at z > 1. 
For example, the hyperluminous galaxy IRAS F10214+4724 (L; r ~ 10 13 L & ; Rowan-Robinson et 
al. 1991) at z = 2.32, which is at least as luminous as the most luminous ULIGs observed locally 
and which is gravitationally lensed in the far infrared by factors as high as 30 (Eisenhardt et 
al. 1996) only just made it into the IRAS 0.2 Jy sample (Oliver et al. 1996). Until more sensitive 
far-infrared/submillimeter space telescopes are launched (e.g. SIRTF and FIRST), and prior to 
the successful operation of sensitive ground-based submillimeter detectors, the most productive 
method of searching for high-z ULIGs (if they exist) seemed to be from trying to identify heavily 
reddened sources in optical, and near-infrared color and redshift surveys. Such methods would 



-3- 



rely heavily on comparing the colors of these candidate objects with the UV-optical colors of 
low-z ULIGs. Within just the past year, an exciting new capability has been added with the 
availability of new deep submillimeter observations obtained with the SCUBA camera on the 
James Clerk Maxwell Telescope (JCMT) (e. g. Smail et al. 1997; Barger et al. 1998; Hughes et 
al. 1998). The new deep submillimeter surveys are able to directly reveal which objects are likely 
ULIG candidates, however optical and near-infrared data are still required to obtain redshifts, 
and to determine more precisely the nature of the source of the intense submillimeter emission. 
An important issue will be to determine if the objects found in these surveys by virtue of their 
submillimeter luminosity are similar to local ULIGs. If this is the case, then we can use the 
detailed measurements of local ULIGs as a template for studying the universe at high redshift. 

When we observe high redshift galaxies at optical and near-infrared wavelengths, we are 
studying their rest-frame ultraviolet properties. Therefore, in order to better identify high-z 
ULIGs, we need to know the rest-frame ultraviolet properties of local ULIGs. The ultraviolet 
colors of ULIGs are mostly determined by internal absorption. Absorption can in principle 
produce arbitrarily red galaxies, although the reddest galaxies would also be faint. It is not 
straightforward to model the color of a typical ULIGs as seen at high z. Ground-based optical 
surveys probe rest-frame ultraviolet wavelengths when z > 1. To make color predictions, we would 
have to model the internal extinction of ULIGs at wavelengths where we have observations only 
from inside our own Galaxy. Its reddening properties are probably different from those of typical 
ULIGs (see, e.g., Kinney et al. 1994). This clearly is dangerous. 

The purpose of this paper is to measure the ultraviolet colors (at ~ 2300 and ~ 1400 A) of 
local ULIGs using the Faint Object Camera (FOC) and the Hubble Space Telescope (HST). We will 
combine the FOC measurements of these local ULIGs with ground-based optical and near-infrared 
data to compute the SEDs for these galaxies. We then compute K corrections and hence predict 
the colors of these galaxies as a function of redshift. We can then compare these colors with the 
observed colors of submillimeter-luminous galaxies found in the recent field surveys. 

Section 2 describes our sample and how it was selected. Sections 3 and 4 present the HST 
and ground-based observations and data reduction, respectively. In § 5, we derive the SEDs and 
K corrections. Section 6 discusses the results and the applications to the high-redshift Universe. 
Section 7 is a summary. 

2. THE SAMPLE 

The basic properties of the galaxies we observed are listed in Table 1. To get a meaningful 
average SED given the (probably different) amounts of extinction in different ULIGs, it is desirable 
to have as large a sample as possible. However, ULIGs are faint in the ultraviolet, so we need to 
take deep exposures. As a result, we were able to observe only three objects. 

Our aim was to select an unbiased sample of ULIGs that might best be representative of 
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their high-redshift counterparts. IRAS F12112+0305 and IRAS F22491-1808 were chosen from 
the IRAS Revised Bright Galaxy Sample (Sanders et al. 1999) because they are typical ULIGs 
based on their "cool" 25 - 60 //m colors (/25//6O < 0.3; de Grijp et al. 1985; /2s and /qq are the 
IRAS flux densities in Jy at 25 and 60 /im respectively). Most ULIGs at z ~ (~80%) are 
cool. Rarer "warm" ULIGs like Mrk 231 which have strong evidence for a non-thermal origin for 
the far infrared luminosity were avoided. Our third source VII Zw 031 has been studied in detail 
by Djorgovski et al. (1990), who interpret it as a merger-induced starburst galaxy. Its IRAS 
luminosity is only 1% below the limiting threshold of 10 12 Lq that Sanders et al. (1988a) used 
to define ULIGs; it therefore was not in their ultraluminous sample (also it is at low Galactic 
latitude). However, we selected it here because it is one of the best examples of an IRAS galaxy 
that appears to derive most of its infrared luminosity from an embedded starburst. We did not 
choose Arp 220, the nearest "cool" ULIG, because it is too big to fit in the small (11 arcseconds) 
field of view of the FOC. 

Ground-based R-band images of the ULIGs observed with the FOC are shown in Figure 1. 
In VII Zw 031, there is no morphological evidence at large radii for an interaction, but at small 
radii, the isophotes deviate strongly from elliptical symmetry. The other two galaxies show much 
stronger evidence for a recent merger. In both cases, the images suggest that the progenitors were 
disk galaxies. The high infrared and CO luminosities (Sage & Solomon 1987, Sanders et al. 1991) 
suggest that the progenitors were gas-rich. 

3. HST OBSERVATIONS AND DATA ANALYSIS 

FOC images of all three galaxies were obtained using the F220W filter (A pca k = 2270 A) 
and the F140W filter (A pea k = 1370 A). Dates, exposure times, and measured fluxes (see below) 
are given in Table 2. Two or three exposures were taken of each galaxy with each filter. The 
offsets between them were found to be negligible (images of a single galaxy were always taken 
sequentially), so the images were simply averaged. All images were recalibrated using the 
CALFOC package in STSDAS; the results were identical to those produced by the HST pipeline 
calibration. They show the usual FOC blemishes: vertical stripes on the right side of the images 
due to variations in the camera scan speed and a fingerprint pattern in the geometrically corrected 
image. Typical count rates in our images were 0.001 to 0.01 counts s _1 pixel -1 ; this is sufficiently 
low that we need not worry about the nonlinearity of the FOC detector (see, e.g., Baxter 1993). 
In the rest of this section we describe the images and how we carried out the photometry. 

3.1. VII Zw 031 = IRA S F 50081 + 7936 

Our F220W image of VII Zw031 is shown in Figure 2. The raw image and two deconvolutions 
using different techniques are shown. The most striking feature is the bright crescent of knots 
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in the deconvolved images. The strong similarity of the two deconvolutions argues that the 
morphology is not sensitive to the deconvolution technique. The knots are also visible in the raw 
image but are less pronounced there. They are not seen in optical ground-based images. While 
the knots are interesting - they may be luminous young star clusters - we do not attempt to carry 
out photometry on the deconvolved images. This is because the deconvolutions do not conserve 
flux locally. All photometry is based on the raw images. 

Aperture brightnesses were measured by centering on the center of curvature of the crescent 
of knots; this position is < 1" from the optical center of the galaxy. The growth curve is shown 
in Figure 3. The limiting factor in determining the total brightness is the uncertainty in the 
sky level. Determining the sky brightness directly is not possible because of the FOC blemishes. 
Instead, we measure the sky value at various places well away from the galaxy and plot growth 
curves for the range of resulting sky values. Five such curves are shown in Figure 3. We then 
identify which growth curve converges best to a plateau as far as possible from the adopted center 
of the galaxy while still indide the FOC field of view, and use it to estimate the total counts in 
the galaxy. The other growth curves provide an estimate of the uncertainties. These are large: 
we adopt a total flux of 9400 ± 2000 counts. However, it will turn out that the uncertainties 
due to sky subtraction are small compared to the differences between different ULIGs. Therefore 
sky subtraction uncertainties do not have a major effect on our conclusions about the average 
properties of ULIGs. 

That the growth curves converge to a plateau at all suggests that the part of the galaxy that 
is significantly above the sky at 2300 A is well within our field of view. However, it is possible 
that the structure we see is embedded in a low-surface-brightness component that contributes 
significantly to the total flux but that is never as bright as the rms sky noise. Such a halo would be 
excluded from the above total flux. We estimate as follows how important it could be. We assume 
that this component has an exponential light profile and a scale length of h = 1 kpc. This is a 
typical scale length for a 0.5 L* galaxy (Freeman 1970; L* is the characteristic luminosity of the 
Schechter 1976 luminosity function); 0.5 L* is probably close to the luminosity of the individual 
galaxies that merged to make VII Zw 031 (Djorgovski et al. 1990). We then measure the rms sky 
noise a at various places in the image and then run a detection algorithm (FOCAS, Jarvis & 
Tyson 1981) to look for objects with fluxes at least 3 a above the background and with sizes at 
least as large as the seeing disk (~ 0"1). When we do this and find no objects, this sets a limit on 
the flux within one seeing disk of this exponential component. Combining this limit with the scale 
length then allows us to set a limit on the total flux in this component. The additional luminosity 
that can be hidden in such a component is 1000 counts for the F220W image of VII Zw 031. This 
is much less than the measured counts. 

One complication is that the sky noise varies from place to place in the image. Also, the 
statistics are altered by the fingerprint pattern introduced by the geometric transform. However, 
a in regions with and without the fingerprint differs by only ~ 20 %; the resulting uncertainty of 
± 200 counts is negligible compared to the uncertainty of ± 2000 counts due to the poorly known 
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sky value. Also, a varies from place to place by ~ 20 % even in regions that do not show the 
fingerprint. We therefore just use a value of a that is averaged over a region large enough to have 
substantial contributions from both fingerprint and non-fingerprint areas. 

The calculations above suggest that we cannot hide a significant fraction of the measured 
flux in a low-surface-brightness component that has a scale length of h ~ 1 kpc. However, if this 
component were more extended, the amount of flux that could go undetected would increase 
rapidly with h (it is approximately proportional to h 2 ). A component with h = 3 kpc could 
dominate the total counts. This is not an unreasonably large size; the string of knots that 
we observe has a length of almost twice this. What this means is that if we wish to compute 
total magnitudes, the values we compute directly from the growth curve in Figure 3 might be 
substantially in error. Therefore we shall restrict ourselves to aperture magnitudes here and in the 
rest of this paper. We may still be missing some of the total flux in our apertures if there exists a 
low surface-brightness component, but this is now a far smaller fraction of the total flux, even if 
this component has h = 3 kpc. 

The above paragraph demonstrates that we need to use aperture magnitudes. We now need 
to select a physical size for the aperture in which to measure fluxes. This is chosen to be 6 kpc in 
each galaxy rest frame, which is the minimum size at which all the F220W growth curves have 
converged to a plateau. It is desirable to use the same aperture for all the galaxies, since in the 
subsequent sections of the paper we shall be comparing their SEDs. 

VII Zw 031 was only detected at the 0.45<r level in the F140W images. These images are 
therefore best regarded as just giving flux limits. A similar aperture analysis to that discussed 
above (see Figure 4) gives a total flux of 250 ± 550 counts. The contribution from a possible 
low-surface-brightness component is proportionately much larger here than for the F220W data. 
We estimate that the total counts that can be hidden in an h = 1 kpc component are 3700, more 
than 10 times the measured counts. 

For both images, the counts are then converted to units of flux density fx using the standard 
FOC calibration equations. We assume that the spectrum is flat across the bandpass (this 
assumption will be examined in detail and corrected for in §5.4). The resulting flux densities 
within the 6 kpc aperture are listed in the final column of Table 2. 

3.2. IRAS F12112+0305 

The F220W image of IRAS F12112+0305 is displayed in Figure 5. The F220W-band emission 
is more diffuse than it is in VII Zw 031; there is no evidence here for luminous star clusters. 
The emission comes from two concentrations; neither is centered near the optical centers of the 
interacting galaxies. Their center positions are shown in Figure 5; they are close to the edge of the 
field of view of the detector. Clearly we are more likely to miss flux at large radii than we were 
for VII Zw 031. In Section 5.3, we discuss how such missing flux would affect the calculation of the 
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SED and the K corrections. 

Figure 6 shows growth curves given by aperture photometry like that in §3.1. We centered 
the apertures mid way between the two light concentrations. Adopting the sky value that results 
in a plateau in the growth curves implies a flux of 15, 000 ± 2, 000 counts. 

The F140W images show faint extended emission distributed like the emission in the F220W 
images. A growth curve analysis like the above (see Fig. 7) gives a measured flux of 6000 ± 2000 
counts. 

The counts were then converted to physical flux units; the results are presented in Table 2. 
As for the previous galaxy, these refer to 6 kpc aperture magnitudes. 

3.3. IRAS F22491-1808 

The F220W and F140W images of IRAS F22491-1808 are shown in Figure 8. We detect this 
galaxy at a far higher level of significance than the other two galaxies. The integration time was 
lower by a factor of ~ 10, yet we have approximately the same number of counts. There is some 
evidence for knots, but they contribute a substantially smaller fraction of the 2300 A flux than 
do the knots in VII Zw 031. Figure 8 also shows that the brightest peak in the F220W image is 
offset from the optical center by about 2". We regard it as unlikely that such a large (> 5<r) peak 
at 2300 A would be undetected in our i?-band optical image. We therefore attribute the offset 
to errors in HST pointing (errors this large are not uncommon - see the FOC Data Handbook, 
Baum 1994). There is further evidence for this interpretation: if we overlay the F220W peak on 
the optical peak, then the faint emission that we see in the F220W image to the east of the main 
galaxy exactly coincides with the eastern extension in the optical image (see Figure 1). 

The photometry in this case is complicated because the galaxy falls close enough to the 
west edge of the image that variations in the camera scan speed cause substantial variations in 
the background in the F140W image (see Figure 8b). However, these variations have a definite 
direction (up-down in the image). We therefore perform photometry not using circular apertures 
but using rectangular apertures, where we vary the aperture dimension in the direction orthogonal 
to the fluctuations. When we subtract the sky, these fluctuations will affect the subtraction in 
a systematic way, so that the mean sky values we subtract which then produce a plateau in the 
growth curve will be slightly different from the mean sky values we measure in the main part of 
the image. These effects are not large (< 10%), and they are easily corrected, as above. We then 
find that the growth curves (Figure 9) converge to a plateau; this also suggests that the measured 
flux above sky drops to zero around the edges of the image. The total flux is 9300 ± 1000 counts. 
The F140W image shows a single brightness peak at the same place as the brightest knot in the 
F220W image and a very faint diffuse component that also follows the emission in the F220W 
image. Aperture photometry gives a total flux of 1500 ± 700 counts. 
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Table 2 again lists the fluxes in physical units. 

4. GROUND-BASED OBSERVATIONS AND DATA ANALYSIS 

To compute SEDs and K corrections as input for the identification of ULIGs in deep surveys, 
we need not only ultraviolet measurements but also optical and near-infrared photometry. We 
obtained this mostly from ground-based observations summarized in Table 3. We also used some 
images taken by other observers and photometry from the literature. All of our observations were 
taken at the University of Hawaii 2.2 m telescope. The detector was a thinned Tektronix 2048 x 
2048 CCD for the optical images and the QUIRC 1024 x 1024 HgCdTe array for the near-infrared 
images. The field of view is 7'5 x 7'5 for the CCD observations at the f/10 focus, 3^2x3md2 for 
the QUIRC observations at f/10, and 1'OxlmdO for the QUIRC observations at f/31. These fields 
were easily large enough to contain all of the galaxy light. We used standard data reduction 
techniques for our images; twilight sky flats were used for the optical images and dome flats for 
the near-infrared images. Most of the optical images were taken as three short exposures; at 
least six were taken for the QUIRC images. Cosmic rays and bad pixels in the vicinity of the 
galaxy were then identified by their appearance in only one of the individual images; they were 
removed by hand. Instrumental magnitudes were computed from observations of 20 - 30 standard 
stars, and the photometry was converted to the Johnson UBV - Cousins RI magnitude system of 
Landolt (1992). For the K' observations, the typical number of standards observed per night was 
6, and we used the color conversions of Wainscoat Sz Cowie (1992) to convert to the magnitude 
system of Elias et al. (1982). The zero points are accurate to 1.5% in the optical and 3% in the 
near-infrared. 

We then measured the magnitudes of the galaxies within an aperture that corresponds to 
6 kpc in the rest frame of each galaxy, as outlined in the previous section. This alignment was 
performed as follows. For VII Zw 031 we centered the optical galaxy at the center of curvature of 
the knots - this is also the position we would choose if we were to assume that the HST pointing 
was completely accurate. For IRAS F22491 — 1808, we aligned the images by requiring the UV 
and optical peaks to coincide, as outlined in Section 3.3. For IRAS F12112+0305, the aperture 
was centered mid way between the UV peaks (see Section 3.2), which is approximately half way 
between the optical peaks (see Figure 5). Adopting these centers, we found that the fraction of 
the light above the 3a isophote falling outside the aperture is about 20% or less. This is small 
enough that highly accurate alignment between the optical and UV images is not essential. For 
these bright galaxies, the main source of error in the magnitudes is the then uncertainty in the 
zero point. The aperture magnitudes are presented in Table 3. 

For IRAS F12112+0305 and IRAS F22491-1808, we obtained several exposures in a number 
of passbands (Table 3). Our B-band image of IRAS F12112+0305 taken in 1995 May revealed 
a bright, point-like object in the southern tidal feature that did not appear in any of the other 
images. This is probably a supernova (Trentham 1997). This fell close to the edge of our aperture 
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(it was 7"7 = 10 kpc from the galaxy center). We measured the total flux in the supernova using 
DAOPHOT (Stetson 1987) and subtracted the part of it that fell within the detection aperture 
to derive the magnitude presented in Table 3. The supernova had a total B magnitude of 18.83 
corresponding to an absolute magnitude of Mb = —18.55 if it is in the galaxy. This supports the 
interpretation as a supernova and suggests that it was observed close to its maximum luminosity 
(Trentham 1997). 

For VII Zw 031, we had no photometric images of our own, so we measured the total 
magnitude on a if-band image from Mazzarella et al. (1999). Optical measurements in the gri 
filter system (Thuan h Gunn 1976; Wade et al. 1979) were taken from Djorgovski et al. (1990). 



SPECTRAL ENERGY DISTRIBUTIONS AND K CORRECTIONS 



In this section we use the results of the previous two sections to compute SEDs for the three 
galaxies. We then compute K corrections; these are defined as follows. 

Suppose we observe the light from a galaxy of redshift z and absolute magnitude Mx in some 
filter X. Its apparent magnitude mx(z) is 



m x (z) = M x + n{z) + K x (z). 



(1) 



Here fi(z) is the distance modulus, and Kx(z) is the K correction. It corrects for two effects: (i) 
the redshifted spectrum is stretched through the bandwidth of the filter, and (ii) the rest-frame 
galaxy light that we see through the filter comes from a bluer part of the SED because of the 
redshift. The second effect is normally the larger, depending on the SED. If K(z) is known for 
filters X and Y, then we can predict the colors as a function of redshift: 



mx(z) - m Y (z) = M x - M Y + K x (z) - K Y (z). 



(2) 



Candidate ULIGs can then be identified in photometric surveys on the basis of color information. 



The K correction for filter X is 



Kx{z) = 2.51og 10 



(1 + *) 



/o 00 /A(V)Tx(AQdy 
f °° fx(£- z )T x (\>)d\> 



(3) 



(e.g., Sandage 1995). Here f\ is the flux density of the source, and ?x(A) is the sensitivity of the 
detector plus filter X. In this work, we measure f\. We can then calculate K corrections for these 
galaxies as functions of z for various filters X. A common simplification that is used when the data 
are poorly sampled (e. g., if only broadband colors are available, as in this paper) is to assume 
that the filter X is narrow; then Tx(A) is a heavily peaked function of wavelength A centered on 
some reference wavelength Aq. That is, ?x(A) pa S(X — Aq). Then equation (3) simplifies to 



K x (z) ^ K Xo (z) = 2.5log 



10 



/a(Aq) 
/a(&) 



(4) 
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Equation (4) is also valid if the SED is flat across the bandpass. We shall use it; the optical filters 
for which we compute K corrections are narrow enough that the uncertainties introduced by 
using equation (4) instead of equation (3) are much smaller than the errors inherent in computing 
fx from the HST observations, where the filters are not narrow. In § 5.4, we show how the 
calculation of fx from broad-band colors (as in § 3) requires us to assume a functional form for 
fx within the bandpasses, and we show how this functional form affects the SED. (In §3, we 
assumed a flat spectrum, a = 0, where fx oc X a .) We then show how the functional form affects 
the computation of K corrections using equation (4). We find that assuming a functional form 
results in compensating systematic errors when the SEDs are calculated and then when the K 
corrections are calculated from the SEDs. We also show in § 5.4 that the magnitude of the errors 
made in using equation (4) are negligible compared to the differences between ULIGs. They are 
therefore not important when we wish to decide whether the objects that we see in color surveys 
are ULIGs or not. 

In the rest of this section, we examine the assumptions and corrections we need to make to 
compute SEDs and K corrections. The results are presented and discussed in the next section. 



5.1. Distance Scale 

We assume a standard Friedmann cosmology with a Hubble constant of Hq = 75 km s _1 
Mpc -1 and a cosmological density parameter of VLq = 1. Then fluxes are related to luminosities as 

L x = ^d 2 L (z)f x . (5) 

Here Lx is the luminosity density (the bolometric luminosity is L = f^° L^dA), and di(z) is the 
luminosity distance. For = 1, 

dL(*) = #(! + *-(! + *)'). (6) 



5.2. Galactic Extinction 

The fluxes and magnitudes quoted in §§ 3 and 4 are the observed values; they need to 
be corrected for Galactic extinction. We adopt S-band Galactic extinctions from Burstein & 
Heiles (1982); they are listed in Table 1. To find the extinctions in other optical bandpasses, 
we use the color conversions of Cardelli et al. (1989). At ultraviolet wavelengths, we used the 
extinction curves given by Mathis (1990). For VII Zw 031, the Galactic extinction is quite large 
(Djorgovski et al. 1990). For Ab = 0.36, the extinction in the F220W filter is A ^ « 0.84. The 
extinction at this wavelength is particularly large because of the 2200 A graphite bump (Mathis 
1990). We apply the Galactic absorption corrections to the observed fluxes or magnitudes before 
computing the luminosity using equation (5). Note that these corrections are much smaller than 
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the corrections we would have to make for internal reddening (see Djorgovski et al. 1990) if we 
wanted to determine the intrinsic luminosities of the galaxies at optical wavelengths. Since our 
aim is to compare local and high-z ULIGs, we do not attempt to correct for internal extinction. 

5.3. Spatial Distribution of Flux 

The FOC field of view is only ll'/3xll'/3. This corresponds to 10.9 kpc for VIIZw031, 14.0 
kpc for IRAS F12112+0305, and 14.9 kpc for IRAS F22491-1808. Most of the light of normal 
galaxies is contained within these diameters, but it is not clear that the same is true for ULIGs. 
Large amounts of internal extinction may reduce the fraction of the ultraviolet light that is near 
the galaxy center. This was in essence our motivation for using aperture magnitudes. 

For all of our images, the growth curves converged on constant total brightnesses well within 
the FOC field of view. While this is encouraging, the caveat is that we determined sky values 
within the same images, so we might have subtracted a low-surface-brightness component that 
extends beyond the edges of the FOC images. Ultraviolet observations with a larger field of view 
would be required to measure this light or show that it is unimportant, meanwhile, we caution 
the reader that the ultraviolet parts of the SEDs may be too faint and therefore that the true K 
corrections may be smaller than the ones we compute at wavelengths A and redshifts z such that 
(1 + z)~ 1 \ < 4400 A. Longward of 4400 A, the K corrections are determined only by the optical 
measurements. 

One other complication is that the IRAS beam (2' x 4') is much larger than the field of view 
of the FOC - therefore when we place far-infrared and UV detections on the same SED, we need 
to be aware of the fact that we might be probing spatially distinct regions of the galaxy. This 
problem only disappears if the far-infrared emission from ULIGs comes from a very compact region 
(which is highly likely; see e.g. Lonsdale et al. 1994 for Arp 220) that is the same region from 
which the UV flux comes. In the rest of this work, we only consider parts of the SED blueward of 
the rest-frame ivT-band, so we can safely ignore this effect. 

5.4. The Slope of the SED Within the Bandpass 

Each magnitude measurement is an integral of the flux density over some filter bandpass, 
J* °° /^Tx(A)dA. To get the SED, we need the flux density fx at some reference wavelength Ao; 
we choose this to be the peak transmission wavelength of the filter-telescope-detector system. 
Therefore, we need to adopt some functional form of f\ within the bandpass. We assume a power 
law f\ oc X a . The problem is that a is unknown. Crude estimates can be obtained by using 
the adjacent points to measure the approximate slopes of the SED on either side of the point 
in question. However, these estimates are not necessarily useful, (i) because the blueward and 
redward point often give very different slopes, and (ii) because the distance between points in 
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wavelength space is much bigger that the bandwidth of the filters. A ore quantitative analysis 
is suggested, as follows. Suppose that we assume some value of a = a\ when we compute our 
SEDs using the total measured flux in filter X with reference wavelength Ao,x an d that we do this 
incorrectly because the true value is a<i- Our SED point will then be incorrect by a factor , 
where 

Ix{a) = / -r x (A)dA • (7) 

Now suppose that we compute the K correction for a filter Y using the incorrect SED that we 
just derived. Filter Y has a redder Ao than filter X because of the redshift z: (1 + z) = The 
denominator of the term inside the brackets in equation (3) will now be incorrect by a factor of 
/ = Jx[q2) jy(tti) ' ano - ^ ne K correction will be incorrect by a factor of 2.5 log 10 /. (We assume 
that the numerator of the term inside the brackets in equation (3) has been computed correctly.) 
Therefore, if the transmission curves for X and Y have the same shape, then / = 1 and the K 
corrections will be independant of a. This is almost true for the optical filters. We therefore 
set a = when we compute the SEDs for the optical and K-band filters. For the HST filters, 
complications are introduced by the broad wings of the F220W and F140W transmission curves 
from Ao toward red wavelengths. Our approach is to assume a = here, too, when we compute 
the SEDs. This value of a is suggested by the IRAS F12112+0305 and IRAS F22491-1808 SED 
slopes between the F220W and F140W point. Doing this also has the advantage that the SED is 
determined consistently for our optical and ultraviolet points. 

However, we still need a quantitative measure of how the assumed a affects our results. 
We therefore present is Table 4 the function 1(a) for the HST filters and for a typical optical 
filter. We can then calculate the / for any particular a and thereby determine how much the K 
correction is affected. For the F140W filter, these numbers can be large. The reason is the broad 
peak in the F140W transmission curve redward of the reference wavelength. For example, at 6000 
A, the F140W filter lets in twice as much light as the F220W filter. However, the SEDs suggest 
that a ~ for this point. For the F220W filter, a is less likely to be 0, because the SED in all 
cases drops strongly between the bluest optical measurement and the F220W data point (although 
some of this is due to the 4000 A break). However, 1(a) is much closer to 1 for the F220W filter 
because it's transmission curve is narrower. We investigated a number of possibilites; typical 
errors introduced by our assumption of a = given plausible values of a range from zero to a few 
percent. The largest is for K corrections which use the F220W point of VII Zw 031. Here, the true 
a may be as large as 4, based on the relative positions of the F220W point and the (/-band point 
(this is the bluest optical point in the SED). This would mean that / is about 16 %. While this is 
significant, it is smaller than the uncertainty in the measured counts and the differences in SEDs 
between the galaxies. Therefore the K corrections are accurate enough to permit identification of 
ULIGs in photometric surveys. 

A potentially more serious worry is the contamination of the broadband colors by emission 
lines. However, Figure 6 of Djorgovski et al. (1990) suggests that the optical broadband colors of 
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VII Zw 031 are not dominated by contributions from the lines. Our FOC niters were chosen to be 
free of known emission lines. We therefore do not expect emission lines to be a problem. 

The SEDs are presented in Figures 11 - 13 and in Table 5. The plots show vL v = \L\ as 
a function of frequency u; they are computed assuming a = 0. The 2270 A and 1370 A points 
derived from the HST photometry are shown as filled circles. For VII Zw 031, the 1370 A point is 
a limit (i. e., a 0.45 a detection). 

From the SEDs, we calculate K corrections Kx(z) for various filters X with reference 
wavelength A , where A (l + z)~ l > 1370 A for IRAS F12112+0305 and IRAS F22491-1808, and 
where Ao(l + z)~ l > 2270 A for VII Zw 031. These K corrections are presented in Table 6. 

6. DISCUSSION 

In this section, we first compare the SEDs of the three sample galaxies and we suggest 
a physical interpretation. We then discuss the implication of our results for ground-based 
submillimeter, near-infrared, and optical color and redshift surveys. 

6.1. The SEDs of the Sample Galaxies 

The SEDs of all three galaxies drop between optical and ultraviolet wavelengths, but by 
substantially different amounts (see Figure 14). VII Zw 031 has the steepest drop; the SED 
drops be two orders of magnitude between the g and F220W bands. IRAS F22491— 1808 has 
the shallowest drop, i.e., by less than a factor of 3 between the B and F220W bands. The 
straightforward interpretation is that the three galaxies have different amounts of internal 
extinction from dust. At optical wavelengths, the internal extinction in VII Zw 031 is Eb~v — 0.9 
(Djorgovski et al. 1990). This suggests that the internal extinction at far-ultraviolet (A < 3000A) 
wavelengths should be huge (factors of hundreds are plausible given typical extinction laws). The 
extremely low ultraviolet flux that we observe is therefore not surprising. The ultraviolet light 
from the other galaxies is presumably less absorbed. It is interesting that these three galaxies 
should show such different properties, given that they were selected in an unbiased way so as to 
be typical of their class. 

Note also that VII Zw 031, which has the most steeply dropping ultraviolet SED, has a 
substantial population of candidate young star clusters in the F220W image. Their ultraviolet 
luminosities are ss 10 s L , comparable to or brighter than the luminosities of the super star 
clusters seen in other merging and starburst galaxies (e.g., Whitmore et al. 1993; Conti & Vacca 
1994; Shaya et al. 1994; Whitmore & Schweizer 1995; O'Connell et al. 1995; Meurer et al. 1995; 
Schweizer 1996; Schweizer et al. 1996; Tacconi-Garman, Sternberg, & Eckart 1997). The clusters 
are several kpc from the galaxy center; we presumably see them because the optical depth is 
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low there. But they contribute much of the ultraviolet light. Therefore ULIGs with the internal 
extinction properties of VII Zw 031 but with no such clusters will have SEDs that drop even more 
steeply in the ultraviolet than does the SED of VII Zw 031. Seen at high redshifts, they would be 
even redder than VII Zw 031. 

6.2. Application to High-Redshift 

In the past three years there have been a number of important breakthroughs in the detection 
and identification of galaxies with redshifts 1 < z < 4. Significant new populations include Lyman 
break galaxies (the ultraviolet dropouts of Steidel et al. 1996) and submillimeter-luminous galaxies 
(Smail et al. 1997, Barger et al. 1998, Hughes et al. 1998, Eales et al. 1999). In addition, a few other 
high redshift galaxies like F10214+4724 (Rowan-Robinson et al. 1991) and HR10 (Hu & Ridgway 
1994) have been identified by virtue of their infrared properties. The submillimeter-luminous 
galaxies appear to have the same submillimeter and far-infrared fluxes as would local ULIGs 
(Sanders & Mirabel 1996) if seen at high redshift (Barger et al. 1998). 

It therefore is productive to compare the observed optical colors of the submillimeter-luminous 
galaxies with our predicted optical colors of local ULIGs were they to be seen at high redshift. 
This comparison will help in (i) establishing the connection between local ULIGs and the 
submillimeter-luminous galaxies, (ii) identifying candidate ULIGs at high redshift in color surveys, 
and (hi) investigating possible links between the Lyman break galaxies of Steidel et al. (1996) and 
local ULIGs (and by implication between the Lyman-break galaxies and submillimeter-luminous 
galaxies) . 

The B — R and B — K colors of our sample galaxies placed at different redshifts are 
presented in Figure 15. Elliptical galaxies are shown for comparison purposes. Color-color and 
color-magnitude diagrams are shown in Figure 16. Kinks in the curves are artefacts of poor 
sampling. The general shapes and positions of the curves are accurate, because the errors in our 
SEDs are small. The forms of the color-color curves are easily understood in terms of the shapes 
of the SEDs. For example, consider the curve for IRAS F121 12+0305 in the upper panel of Figure 
16. Initially, the gradient of the curve is shallow: the S-band fluxes drop rapidly with increasing 
z due to the drop of the SED between optical and ultraviolet wavelengths, but the R and K 
measurements are still probing a flat part of the SED and so have similar K-corrections. When 
the i?-band measurements begin to probe wavelengths corresponding to the optical-UV drop, the 
B — R color remains approximately constant, but the B — K color continues to redden. Therefore 
the curves turn upward in the color-color diagram. Eventually, the S-band measurements begin 
to probe a flat part of the SED (between the F220W and F140W measurements), while the 
i?-band measurements continue to probe wavelengths in the region of the optical-UV drop. The 
curves then become horizontal: the B — R color becomes bluer, but the B — K color remains 
approximately constant. The IRAS F22491— 1808 curve is somewhat different: the optical-UV 
drop is sufficiently shallow that the bandwidth part of the K correction — the (1 + z) term inside 
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the brackets in equations (3) and (4) — is proportionally more important. 

The figures suggests that only one of our galaxies (VII Zw 031) can appear redder in B — R 
than typical elliptical galaxies at z = 0.8. This is the redshift at which elliptical galaxies are 
reddest - see Coleman et al. (1980). The B — K colors of high-z ellipticals are redder than those 
of all our sample galaxies at all redshifts. It is therefore clear that we cannot hope to obtain 
complete samples of ULIGs at high-z using color surveys alone. For example, we could not have 
found the hyperluminous high-z galaxy IRAS F10214+4724 this way. Nevertheless the reddest 
ULIGs should easily be detectable at z > 1, as they would have colors redder than any normal 
galaxies. For example, very red faint galaxies have been found by Hu Sz Ridgway (1994), one of 
which (HR10) has since been detected by SCUBA (Cimatti et al. 1999). 

A comparison with the V — I colors of the SCUBA sources detected by Smail et al. (1997, 
1998) is given in Figure 17. We use V — I colors as opposed to the B — K colors presented in Figure 
16 because i^-band colors are not yet available for a significant number of SCUBA sources. Note 
that the Smail et al. (1998) galaxies are seen through low redshift clusters and so are magnified 
by as much as a factor of ~ 2 — 5 through gravitational lensing. This means that the we need to 
shift the points in Figure 17 horizontally to the right by up to one magnitude when comparing the 
points with the lines. The V — I colors are only weakly affected by gravitational lensing, if at all. 

We see that VII Zw 031 at high redshift would appear redder in V — I than any of the Smail 
et al. galaxies; on the other hand all the Smail et al. galaxies (except the brightest two, which 
are cD galaxies) have V — I colors close (within one magnitude) to those of IRAS F12112+0305 
and IRAS F22491— 1808 seen at high redshift. The high-redshift Lyman break galaxies, recently 
discovered by Steidel et al. (1996), are also systematically bluer (in R — K) than our sample 
galaxies would appear if seen at the same redshift (see Table 7). We use R — K colors for this 
comparison (as opposed to B — K) because our measurements do not probe far enough into the 
ultraviolet to permit us to predict B magnitudes for our sample galaxies placed at the very high 
redshifts of the sample of Lyman break galaxies used in Table 7. 

One implication of these results is that current studies of large optically-selected samples of 
field galaxies that attempt to determine cosmological star formation rates (e.g. Madau et al. 1996) 
will be missing very red ULIGs like VII Zw 031. Furthermore, if star formation rates are derived 
solely on the basis of rest-frame ultraviolet colors, they may then be incorrectly determined 
(Hughes et al. 1998) for galaxies like IRAS F12112+0305 and IRAS F22491-1808, which have 
optical colors of normal galaxies seen at high redshift. If ULIGs are a significant population at 
high redshift, as they appear to be (Smail et al. 1997, Barger et al. 1998, Hughes et al. 1998, Eales 
et al. 1999), then these optical studies will clearly underestimate the star formation rate of the 
universe at high redshift. ore recently, the need to correct optically-derived star formation rates 
for dust extinction has been recognized (Calzetti 1999, Pettini et al. 1999). Deep submillimeter 
surveys like those currently being carried out with SCUBA will in the future determine how 
necessary such corrections are and whether the corrections mostly come from a small amount of 



-16- 



extinction in average galaxies or from a few extreme galaxies being misidentified or issed from the 
samples altogether. 

7. SUMMARY 

We have presented HST FOC images of three low redshift ULIGs and used the measured 
fluxes to compute SEDs for the galaxies. VII Zw 031, IRAS F12112+0305, IRAS F22491-1808) 
were selected to span the range of properties of ULIGs; they have very different ultraviolet SEDs. 
All have SEDs that drop between optical and ultraviolet wavelengths, but the drops range from a 
factor of - 3 in IRAS F22491-1808 to a factor of - 100 in VII Zw 031. We interpret this as being 
due to different internal extinctions in the galaxies. A high ultraviolet extinction in VII Zw 031 
is also suggested by extrapolating the optical internal extinction measured by Djorgovski et 
al. (1990). Implications of our results are as follows: 

1. Only one of our three galaxies (VII Zw 031) has an SED that drops steeply enough to 
ultraviolet wavelengths that its B — R color as seen at high z is redder than the B — R colors 
of elliptical galaxies seen at any z. The other two galaxies have B — R colors that are normal 
for faint field galaxies. Therefore only the reddest ULIGs can be found at high z on the basis 
of optical and near-infrared colors alone. The implication is then that submillimeter surveys, 
like those with SCUBA, are required to obtain complete samples of ULIGs at high z. 

2. Galaxies like VII Zw 031 seen at high redshift may not appear in any current optically 
selected field galaxy sample because their rest-frame UV fluxes are too attenuated. Galaxies 
like IRAS F12112+0305 and IRAS F22491-1808 would appear in these samples, but their 
star-formation rates, if derived solely from their rest-frame ultraviolet fluxes, may be serious 
underestimates. Therefore, any conclusions about the star formation rate of the Universe 
that are derived from optically selected samples are subject to the caveat that ULIGs may 
be missing or not correctly identified. The submillimeter surveys currently in progress will 
tell us how serious this worry is. 

3. VII Zw 031 seen at high redshift would have a redder V — I colour than any of the SCUBA 
sources studied at optical wavelengths by Smail et al. (1998). However, most of the Smail et 
al. galaxies have similar V — I colors to either IRAS F12112-0305 or IRAS F22491-1808 
seen at z > 1. The Lyman break galaxies in the sample of Pettini et al. (1998) have R — K 
colors that are systematically bluer than those of our sample galaxies if seen at the same 
redshift. 
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FIGURE CAPTIONS 



Fig. 1. — Ground-based R-band images of our sample galaxies. Each field is 51" x 51"; this 
corresponds to 49 x 49 kpc for VII Zw 031, 63 x 63 kpc for IRAS F12112+0305, and 67 x 67 kpc 
for IRAS F22491-1808. All images in this paper have north up and east to the left. The VII Zw 031 
image is from the data of Mazzarella et al. (1999); the other images are from the University of Hawaii 
2.24 m telescope (see §4). 

Fig. 2.— The HST FOC image of VII Zw 031 taken through the F220W filter. The optical center 
of the galaxy is identified with a +. The figure shows a section of the image about 7"3 = 7 kpc on 
a side. The left panel shows the raw data. The other two panels show deconvolved data. For all 
of the deconvolutions, the point-spread functions (PSF) were generated using the HST TinyTim 
software for a star with B — V = 1.1. The center panel shows the data after 20 iterations of Lucy 
(1974) - Richardson (1972) deconvolution. The right panel shows the data after maximum entropy 
deconvolution using the STSDAS MEMSYS-5 software with the MEM front end (N. Weir 1991, 
private communication). It is the mean of two images, which are the raw data deconvolved with 
the PSF described above and using intrinsic correlation functions of correlation length 0"02 and 
0"16 (the PSF has FWHM ~ 0"1). An explanation of the meaning and use of intrinsic correlation 
functions and a discussion of the advantages of the above approach are given in detail in the MEM 
manual. 

Fig. 3. — Growth curves for the F220W image of VII Zw 031 based on aperture photometry carried 
out as described in the text. The sky values for the six growth curves are 2170, 2190, 2210, 2230, 
and 2250 counts arcsec -2 from the top down. The heavy line plateaus over the widest range in 
radius; we adopt it and use the others to estimate the errors. The resulting total brightness of 
VII Zw 031 at 2300 A is 9400 ± 2000 counts. 

Fig. 4. — Growth curves for the F140W image of VII Zw 031. The corresponding sky values are 
5170, 5210, 5250, 5290, and 5330 counts arcsec -2 from the top down. 

Fig. 5.— The F220W image of IRAS F12112+0305. The box is a square of side 9"l = 12.1 kpc. 
The two crosses mark the centers of the two galaxies that make up this interacting system. The 
fingerprint pattern introduced by the geometric distortion correction is visible at a low level in this 
image. 

Fig. 6.— Growth curves for the F220W image of IRAS F12112+0305. The corresponding sky 
values are 1550, 1570, 1590, 1610, 1630, and 1650 counts arcsec -2 from the top down. 

Fig. 7. — Growth curves for the F140W image of IRAS F12112+0305. The corresponding sky 
values are 4545, 4590, 4630, 4650, 4670, 4710, and 4755 counts arcsec -2 from the top down. 
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Fig. 8.— (a) The F220W image of IRAS F22491-1808. The box is a square of side 11.3 arcseconds 
= 16.0 kpc. The + marks the optical center of the galaxy; (b) the F140W image of IRAS 
F22491— 1808. The scale and orientation are the same as in (a). 

Fig. 9. — Growth curves for the F220W image of IRAS F22491— 1808. Rectangular apertures were 
used as described in the text; r represents the distance from the center to the edge of the aperture 
in the direction orthogonal to that for which large-scale fluctuations in the detector are seen. The 
width of the detector in that direction is then (2r + 1) arcseconds. The width of the aperture in 
the orthogonal direction (i. e., that of the fluctuations) is fixed at 8''9. The sky values are 166, 176, 
186, 196, and 206 counts arcsec -2 from the top down. 

Fig. 10. — Growth curves for the F140W image of IRAS F22491-1808. Radius has the same 
meaning as in Figure 9. The sky values are 920, 940, and 960 counts arcsec -2 from the top down. 

Fig. 11. — The rest-frame spectral energy distribution of VII Zw 031. The open symbols represent 
the IRAS 12, 25, 60, and 100 (im measurements (the uncertainties in the IRAS measurements are 
smaller than the symbols) and the ground-based g, r, i, and K-band measurements. The filled 
symbols and error bars are the HST FOC results. 

Fig. 12. — The rest-frame spectral energy distribution of IRAS F12112+0305. Open and closed 
symbols have the same meaning as in Figure 11. The optical and near- infrared bandpasses are B, 
V, R, /, and K' . 

Fig. 13. — The rest-frame spectral energy distribution of IRAS F22491-1808. Open and closed 
symbols have the same meaning as in Figure 11. The optical and near-infrared bandpasses are B, 
V, R, I, and K. 

Fig. 14. — The SEDs of our sample galaxies, all (arbitrarily) normalized to have the same luminosity 
as F10214+4724 at x 10 15 Hz. The symbols have the same meanings as In Figures 11 to 13. The 
error bars for all the points are given in Figures 11 to 13; note that the bluemost point for VII Zw 031 
represents only a 0.45<r detection. 

Fig. 15. — The B — R and B — K colors of our sample galaxies as a function of redshift. These 
were computed using the K corrections in Table 6. The colors of Mb = — 20 elliptical galaxies are 
taken from the measurements of Coleman et al. (1980). We use the near-infrared K corrections for 
elliptical galaxies of Huang et al. (1997). 

Fig. 16. — The upper panel is a color-color diagram showing how our sample galaxies would appear 
at different redshifts. The elliptical galaxy is assumed to have Mb = —20 at z = 0. Each of the 
curves is labelled with two numbers. The end labelled "0" represents the galaxy at z = 0. We 
present each curve to the highest redshift that is constrained by our SEDs; this redshift is the label 
at the other end of the curve. We use F140W data for IRAS F12112+0305 and IRAS F22491-1808 
only. The lower panel is a corresponding color-magnitude diagram. The curves and symbols have 
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the same meaning as in the upper panel. Elliptical galaxies are not included on this panel; they 
can occupy almost any region of this plot, depending on their absolute magnitude and redshift. In 
both panels, IRAS F10214+4724 (Rowan-Robinson et al. 1993) is indicated as a cross. 

Fig. 17. — The V — I versus / color-magnitude diagram for our sample galaxies seen at different 
redshifts. The curves have the same meanings as in the lower panel of Figure 16. The filled 
circles represent the colors and magnitudes for the sample of SCUBA sources studied by Smail et 
al. (1998). 
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TABLE 1 

Galaxy Properties 1 



Name 


RA. Dec. 


RA. Dec. 








z 


A B 








12/im 


25/im 


60/im 


100/mi 








B1950.0 


J2000.0 


Jy 


Jy 


Jy 


Jy 




mag 


VII Zw 031 


05:08:07.6 79:36:44 


05:16:38.9 79:40:11 


0.21 


0.58 


5.58 


9.62 


0.0542 


0.36 


IRAS F12112+0305 


12:11:13.7 03:05:15 


12:13:47.3 02:48:34 


0.12 


0.51 


8.50 


9.98 


0.0723 


0.04 


IRAS F22491-1808 


22:49:08.7 -18:08:24 


22:51:49.0 -17:52:27 


0.05 


0.55 


5.44 


4.45 


0.0777 


0.03 



Positions from the IRAS FSC (Moshir et al. 1992). Redshifts and fluxes from the RBGS (Sanders et al. 1999) 
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TABLE 2 

HST FOC Ultraviolet Observations 



Galaxy 


Filter 


Date 


Exposure 


(fx)obs 










sec 


erg cm~ 2 s _1 A~ 


i 


VII Zw 031 


F220W 


8 Nov. 1993 


5390 


(1.28 ±0.23) x 10 


-16 




F140W 


8 Nov. 1993 


8090 


(1.31 ±2.88) x 10 


-17 


IRAS F12112+0305 


F220W 


11 Jun. 1993 


2393 


(2.52 ±0.34) x 10 


-16 




F140W 


11 Jun. 1993 


4794 


(3.54 ± 1.18) x 10" 


-16 


IRAS F22491-1808 


F220W 


8 May 1993 


353.4 


(1.06 ±0.12) x 10" 


-15 




F140W 


8 May 1993 


593.4 


(7.14 ±3.25) x 10" 


-16 
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TABLE 3 

Optical and Near-Infrared Observations 



Galaxy 


Filter 


Run 1 


Exposure 
min 


Magnitude 


IRAS F12112+0305 


B 


2 


8 


17.66 




V 


4 


5 


16.94 




R 


4 


15 


16.40 




I 


4 


5 


15.63 




K' 


1 


6 


12.88 


IRAS F22491-1808 


B 


4 


10 


17.08 




V 


4 


10 


16.42 




R 


4 


20 


16.02 




I 


4 


5 


15.55 




K 


3 


20 


13.01 



Observing runs: 1. 11-13, Mar 1995 (UH 2.2 m with QUIRC 1024 x 1024 HgCdTe array at f/10); 2. 30 Apr-3 
May 1995 (UH 2.2 m with Tektronix 2048 x 2048 CCD at f/10); 3. 14-16 Jul 1995 (UH 2.2 m with QUIRC at f/31); 
4. 1 Aug 1995 (UH 2.2 m with Tektronix CCD at f/10). 
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TABLE 4 

1(a) for F220W, F140W, and R-band Filters 1 



a 


F220W 


F140W 

7(a) 


R-band 


—4 o 


1 01 

J. .U J. 


68 


Q7 


-3.5 


1.00 


0.70 


0.97 


-3.0 


0.99 


0.72 


0.97 


-2.5 


0.99 


0.75 


0.98 


-2.0 


0.98 


0.78 


0.98 


-1.5 


0.98 


0.82 


0.99 


-1.0 


0.99 


0.87 


0.99 


-0.5 


0.99 


0.93 


0.99 


0.0 


1.00 


1.00 


1.00 


0.5 


1.01 


1.09 


1.01 


1.0 


1.03 


1.22 


1.01 


1.5 


1.05 


1.39 


1.02 


2.0 


1.07 


1.64 


1.03 


2.5 


1.10 


2.01 


1.03 


3.0 


1.14 


2.57 


1.04 


3.5 


1.18 


3.46 


1.05 


4.0 


1.23 


4.91 


1.06 


4.5 


1.29 


7.33 


1.07 


5.0 


1.37 


11.45 


1.08 



J The calculation for F220W and F140W filters assumes 
that the detector is the FOC on the HST. The calculation 
for the R-band filter assumes a Mould filter (A = 6400 A) 
using a thinned Tektronix CCD on the UH 2.24 m tele- 
scope. 
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TABLE 5 

Spectral Energy Distributions 



Galaxy 




^rest frame 

Hz 


(^^tOrest frame 

erg s _1 




Source 


VTT Tttt HQ1 
Vll ZjW Uol 


i nn 


o if; v i r\12 
O.10 X 1U 


1 00 \y 1 n45 

i.yz x iu 




IRAS 




fin 


K 07 ss i ni2 
O.Z/ X 1U 


1 SI ss 1 n45 
l.ol x iu 




IRAS 




O £ 


i o£ \/ i ni3 
1.Z0 X 1U 


/j co v i n44 
4.0o X IU 




IRAS 




1 O 
1Z 


O £Q \/ i ni3 
Z.Oo X 1U 


o oo w i n44 
o.Zo X IU 




IRAS 




z.z 


1 A A \y 1 Ol4 

1.44 X 1U 


1 OQ \y 1 n44 
l.Zo X IU 


Sanders 1996b 




U.oZ 


o.oO X 1U 


k t^n \s i n43 

O.OU X IU 


Djorj 


?ovski et al. 1990 




U.00 


a on w i f)14 
4.oo X 1U 


o no \s i n43 
y.UZ X IU 


Djorj 


jovski et al. 1990 




U.4y 


0.4Z X 1U 


1 OS ss i n44 
l.Zo X IU 


Djorj 


jovski et al. 1990 




0.23 


1.39 x 10 15 


(2.87 ± 0.06) x 10 42 




HST FOC 




U.14 


o qi v 1 nis 
Z.ol X 1U 


(n no -i- A r~i \ \/ i n4i 
^z.uy±4.oij x iu 




HST FOC 


IxlAo r IZllZ+UoUO 


i nn 


i oo \/ i ni2 
o.ZZ X 1U 


QO ss i n45 
Z.oZ X IU 




IRAS 




^n 
ou 


0.00 X 1U 


/I 71 v 1 n45 
4. / 1 X IU 




11 LjL u 




O £ 

ZO 


i oo \s i ni3 
i.zy x iu 


oo \y i n44 

o.yy x iu 




IRAS 




1 
1Z 


o (\q \y i ni3 

Z.Oo X 1U 


Q 00 \y 1 n44 

o.yz x iu 




IRAS 




Z.Z 


i t^o \y i ni4 
l.OZ X iu 


k (\r\ \y i n43 

O.OU X IU 




UH 2.24 m 




n o 
u.y 


Q E^7 \/ i ni4 
O.O/ X IU 


A 1 A \y 1 n43 
4.14 X IU 




UH 2.24 m 




n 7 


/i co v i ni4 
4.oy x iu 


o 07 w i n43 
o.oi X IU 




UH 2.24 m 




n 

U.00 


k co v i ni4 
O.oy X IU 


o oc v 1 n43 
o.oO X IU 




UH 2.24 m 




U.44 


7 qi Y i n!4 
/.ol X IU 


o 7C v 1 n 4 3 
Z. / X IU 




UH 2.24 m 




0.23 


1.41 x 10 15 


(6.07 ±0.81) x 10 42 




HST FOC 




0.14 


2.35 x 10 15 


(5.13 ± 1.71) x 10 42 




HST FOC 


IRAS F22491-1808 


100 


3.23 x 10 12 


1.79 x 10 45 




IRAS 




60 


5.39 x 10 12 


3.44 x 10 45 




IRAS 




25 


1.29 x 10 13 


8.74 x 10 44 




IRAS 




12 


2.69 x 10 13 


2.93 x 10 44 




IRAS 




2.2 


1.47 x 10 14 


6.01 x 10 43 




UH 2.24 m 




0.9 


3.59 x 10 14 


5.01 x 10 43 




UH 2.24 m 




0.7 


4.62 x 10 14 


5.34 x 10 43 




UH 2.24 m 




0.55 


5.88 x 10 14 


5.96 x 10 43 




UH 2.24 m 




0.44 


7.35 x 10 14 


3.87 x 10 43 




UH 2.24 m 




0.23 


1.42 x 10 15 


(2.88 ±0.31) x 10 43 




HST FOC 




0.14 


2.36 x 10 15 


(1.17 ± 0.55) x 10 43 




HST FOC 
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TABLE 6 

K-CORRECTIONS 



VII Zw 031 IRAS F12112+0305 IRAS F22491-1808 



z 






















J^-B 




J^-R 






0.1 


0.51 


-0.13 


-0.23 


0.09 


0.09 


0.13 


0.03 


0.06 


0.05 


0.03 


0.17 


0.01 


-0.03 


0.01 


0.02 


0.2 


0.97 


-0.11 


-0.37 


0.05 


0.17 


0.35 


0.11 


0.07 


0.12 


0.07 


0.21 


0.19 


-0.07 


-0.01 


0.04 


0.3 


1.40 


0.31 


-0.48 


-0.14 


0.25 


0.54 


0.18 


0.07 


0.20 


0.10 


0.25 


0.36 


-0.11 


-0.03 


0.06 


0.4 


1.80 


0.71 


-0.58 


-0.32 


0.32 


0.73 


0.32 


0.09 


0.25 


0.13 


0.28 


0.45 


-0.10 


-0.06 


0.07 


0.5 


2.16 


1.07 


-0.66 


-0.46 


0.38 


0.90 


0.49 


0.16 


0.25 


0.16 


0.32 


0.48 


0.05 


-0.09 


0.09 


0.6 


2.51 


1.42 


-0.31 


-0.53 


0.44 


1.06 


0.65 


0.22 


0.25 


0.18 


0.35 


0.52 


0.18 


-0.12 


0.10 


0.7 


2.83 


1.74 


0.01 


-0.61 


0.50 


1.21 


0.80 


0.28 


0.26 


0.21 


0.38 


0.54 


0.31 


-0.15 


0.12 


0.8 


3.14 


2.05 


0.31 


-0.69 


0.55 


1.35 


0.94 


0.41 


0.28 


0.23 


0.41 


0.57 


0.35 


-0.13 


0.13 


0.9 


3.43 


2.34 


0.60 


-0.76 


0.60 


1.48 


1.07 


0.55 


0.33 


0.25 


0.43 


0.60 


0.38 


-0.02 


0.14 


1.0 


3.70 


2.61 


0.88 


-0.58 


0.65 


1.61 


1.21 


0.67 


0.38 


0.27 


0.46 


0.62 


0.40 


0.09 


0.15 


1.1 


3.96 


2.87 


1.14 


-0.32 


0.69 


1.71 


1.33 


0.80 


0.43 


0.29 


0.50 


0.65 


0.43 


0.19 


0.16 


1.2 


4.21 


3.12 


1.38 


-0.07 


0.74 


1.73 


1.44 


0.91 


0.47 


0.30 


0.59 


0.67 


0.45 


0.29 


0.17 


1.3 




3.36 


1.62 


0.17 


0.78 


1.74 


1.55 


1.02 


0.58 


0.32 


0.67 


0.69 


0.47 


0.32 


0.18 


1.4 




3.58 


1.85 


0.39 


0.82 


1.76 


1.66 


1.13 


0.69 


0.33 


0.75 


0.71 


0.49 


0.34 


0.19 


1.5 




3.80 


2.07 


0.61 


0.86 


1.77 


1.76 


1.23 


0.79 


0.35 


0.83 


0.73 


0.51 


0.36 


0.20 


1.6 




4.01 


2.28 


0.82 


0.89 


1.79 


1.85 


1.32 


0.89 


0.39 


0.91 


0.75 


0.53 


0.38 


0.20 


1.7 




4.21 


2.48 


1.02 


0.93 


1.80 


1.87 


1.42 


0.98 


0.42 


0.98 


0.82 


0.55 


0.40 


0.18 


1.8 






2.67 


1.22 


0.85 


1.82 


1.88 


1.51 


1.07 


0.45 


1.05 


0.89 


0.57 


0.41 


0.17 


1.9 






2.86 


1.40 


0.76 


1.83 


1.89 


1.60 


1.16 


0.49 


1.12 


0.96 


0.59 


0.43 


0.17 


2.0 






3.04 


1.59 


0.68 


1.84 


1.91 


1.68 


1.24 


0.52 


1.18 


1.03 


0.60 


0.45 


0.16 


2.1 






3.21 


1.76 


0.60 


1.85 


1.92 


1.76 


1.32 


0.54 


1.24 


1.09 


0.62 


0.46 


0.15 


2.2 






3.38 


1.93 


0.53 


1.86 


1.93 


1.84 


1.40 


0.57 


1.31 


1.15 


0.63 


0.48 


0.14 


2.3 






3.55 


2.09 


0.45 


1.88 


1.94 


1.92 


1.48 


0.58 


1.37 


1.21 


0.65 


0.49 


0.13 


2.4 






3.71 


2.25 


0.39 


1.89 


1.95 


1.93 


1.55 


0.58 


1.43 


1.27 


0.70 


0.51 


0.11 


2.5 






3.86 


2.41 


0.35 


1.90 


1.96 


1.94 


1.63 


0.58 


1.48 


1.32 


0.76 


0.52 


0.10 


2.6 








2.56 


0.31 


1.91 


1.97 


1.95 


1.70 


0.58 


1.54 


1.38 


0.81 


0.54 


0.08 


2.7 








2.70 


0.27 




1.98 


1.96 


1.76 


0.58 




1.43 


0.86 


0.55 


0.07 


2.8 








2.85 


0.24 




1.99 


1.97 


1.83 


0.58 




1.48 


0.91 


0.56 


0.06 


2.9 








2.98 


0.20 




2.00 


1.98 


1.89 


0.58 




1.53 


0.96 


0.58 


0.04 


3.0 








3.12 


0.17 




2.01 


1.99 


1.96 


0.58 




1.58 


1.01 


0.59 


0.03 


3.1 








3.25 


0.14 




2.02 


2.00 


2.02 


0.59 




1.63 


1.06 


0.60 


0.02 


3.2 








3.38 


0.10 




2.03 


2.01 


2.08 


0.61 




1.68 


1.11 


0.61 


0.06 


3.3 








3.51 


0.07 




2.04 


2.02 


2.11 


0.63 




1.72 


1.15 


0.64 


0.11 


3.4 








3.63 


0.04 




2.05 


2.02 


2.12 


0.65 




1.77 


1.20 


0.68 


0.16 


3.5 








3.75 


0.02 




2.05 


2.03 


2.13 


0.67 




1.81 


1.24 


0.72 


0.20 


3.6 










0.13 






2.04 


2.14 


0.69 






1.28 


0.77 


0.25 


3.7 










0.25 






2.05 


2.14 


0.72 






1.33 


0.81 


0.29 


3.8 










0.37 






2.06 


2.15 


0.73 






1.37 


0.85 


0.34 


3.9 










0.48 






2.06 


2.16 


0.75 






1.41 


0.89 


0.38 


4.0 










0.58 






2.07 


2.17 


0.77 






1.44 


0.93 


0.42 


4.1 










0.69 






2.08 


2.18 


0.79 






1.48 


0.97 


0.47 


4.2 










0.79 






2.08 


2.18 


0.84 






1.52 


1.00 


0.49 


4.3 










0.89 






2.09 


2.19 


0.88 






1.56 


1.04 


0.50 


4.4 










0.99 






2.09 


2.20 


0.93 






1.59 


1.08 


0.51 


4.5 










1.09 






2.11 


2.20 


0.97 






1.63 


1.11 


0.52 
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TABLE 7 

Comparison with R-K Colors of Lyman Break Galaxies at High Redshift 



Galaxy 


z 1 


R-K 1,2 


R K 
VII Zw 031 at z 


R K 

IRAS F12112-0305 at z 


R K 

IRAS F22491-1808 at z 


Q0000-263D6 


2.96 


2.28 




4.37 


3.37 


Q0201+113C6 


3.05 


2.36 




4.37 


3.44 


Q0201+113B13 


2.17 


2.79 


6.10 


4.23 


2.91 


B2 0902+343 C6 


3.08 


3.07 




4.37 


3.46 


DSF 2237+116 C2 


3.32 


3.11 




4.35 


3.46 



J Data from Pettini et al. (1998a) 

2 The R filter is similar to the Kron-Cousins R filter used elsewhere in this paper (A e ff = 0.69 /tm; Steidel & Hamilton 
1992) 
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